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Abstract 

We obtain the band edge eigenfunctions and the eigenvalues of solvable periodic potentials using 
the quantum Hamilton - Jacobi formalism. The potentials studied here are the Lame and the asso- 
ciated Lame which belong to the class of elliptic potentials. The formalism requires an assumption 
about the singularity structure of the quantum momentum function p, which satisfies the Riccati 
type quantum Hamilton - Jacobi equation, p 2 — ih-^p = 2m(E — V(x)) in the complex x plane. 
Essential use is made of suitable conformal transformations, which leads to the eigenvalues and 
the eigenfunctions corresponding to the band edges in a simple and straightforward manner. Our 
study reveals interesting features about the singularity structure of p, responsible in yielding the 
band edge eigenfunctions and eigenvalues. 
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I. INTRODUCTION 



As is well - known, periodic potentials play a significant role in condensed matter physics. 
Recently, optical lattices have manifested in Bose - Einstein condensates (BEC)P leading 
to the analysis of the excitation spectra in such systems. The energy spectrum of periodic 
potentials is unique due to the existence of energy bands and the solutions of the Schrodinger 
equation have the Bloch form: 

ij)(x) = u(x) exp(ik.x). (1) 

Here u(x) has the periodicity of the potential function. These features of the periodic 
potentials are usually illustrated using the Kronig - Penney modelP in which the crystal 
lattice is approximated by a periodic array of square wells. Imposing the Bloch condition 
on the solutions of the Kronig- Penney model, one is led to a transcendental equation, which 
has to be solved in order to obtain the band edge solutions. 

Of late, there has been considerable interest in the study^El and construction^-^ of new 
periodic potentials. The periodic lattices in BEC has also renewed interest in this areaP 
The number of exactly solvable (ES) potentials is very limited and in that, exactly solvable 
periodic potentials constitute a very small number. Quite sometime back it was shown by 
Scarf that, V(x) = A/ sin 2 re, usually treated in a finite domain, can lead to band structure 
for -1/4 < A < oE@. In this paper, we study the Lame and the associated Lame potentials 
which belong to the class of elliptic potentials which are usually expressed in terms of the 
Jacobi elliptic functions sn (x, m), cn (x, m) and dn (x, m). The parameter m is known as the 
elliptic modulus whose value lies between and 1. By applying supersymmetric quantum 
mechanics and group theoretical techniques to the family of elliptic potentials, many 
new ES and QES potentials have been constructed ^ HE l 

The Lame potential 

V(x) = j(j + l)msn 2 (x, m), (2) 

with j being an integer, has 2j + 1 bands followed by the continuum. This potential has 
been extensively studieeff^"^. It is interesting to note that the Schrodinger equation with 
the Lame potential 



d 2 ip(x) h 



2 



dx 2 2m 



H (E — j(j + l)m sn (x, m))ijj(x) = 0, (3) 
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is a result of separating the Laplace equation in the ellipsoidal coordinates and is referred to 
as the Lame equatior l^^l This has unique analytical properties^!, which render it useful 
in the study of topics ranging from astrophysics E3-G21 to condensed matter physics^! For 
example one comes across Lame equation in problems related to the early epoch of the 
universe like the distance red-shift relations and quantum vacuum fluctuations.^^ The 
Lame equation occurs in the study of bifurcations in chaotic Hamiltonian systemJ^EHl it 
was shown by Finkel eiaZ^^that the matter modes equation, corresponding to the most 
general inflation potential in the Minkowski space reduces to the Lame form. 
In comparison, the associated Lame potential 

in \ 2/ \ cn 2 (x,m) 

V(x)=pmsn {x,m) + qm — 57 r> (4) 

dn [x, to) 

with p = a(a + 1) and q = b(b + 1), did not have a systematic study till now and has a few 
scattered references in the mathematical literature. This potential is ES when a = b = j , 
with j being an integer and QES when a 7^ b with a, b being real. Recently, a systematic 
study of the band structure of this potential, for various cases of a and b, has been done 
by Khare and Sukhatmd^ and the algebraic aspects of the QES associated Lame potential 
have been investigated irP^H^. 

The band edge eigenf unctions and the eigenvalues of both these potentials have been 
found using the conventional techniques of solving the ordinary differential equation J-^^ i 
Here, we propose an alternate method through the quantum Hamilton - Jacobi (QHJ) 
formalism, where one can obtain the solutions, without solving the differential equations 
explicitly. The QHJ formalism was shown to yield bound state energies and wave functions 
for a large number of ES models. The tools used for our earlier studies are sufficient to 
derive the results for the periodic potentials also. This method requires making a guess 
about the singularity structure of the logarithmic derivative of the wave function ip and using 
elementary tools from complex variables to fix the QMF completely. Once the form of the 
QMF is obtained getting the band edge wave functions and eigenvalues is straightforward. 

In the next section, we give a brief summary of the QHJ formalism and the steps involved 
in obtaining the solutions for any general potential. In sec III, we show how one can obtain 
the form of the wave functions for the Lame potential for a general j, followed by sec IV, 
where we take up a special case of j = 2 and obtain explicit band edge wave functions and 
energies. In sec V, the form of the wave functions for the general associated Lame potential 
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(ES case) is derived. Associated Lame potential, with a = b = 1, is worked out in sec VI, 
followed by the last section containing the concluding remarks. 



II. QUANTUM HAMILTON- JACOBI FORMALISM 

The QHJ formalism, was initiated by Leacock and Padgett P^^ in 1983 - 84. Using 
this formalism on a host of ES models, Bhalla et. were successful in obtaining the 

corresponding energy eigenvalues. The method when applied to the QES models, yielded the 
quasi exactly solvability conditional With a slight modification of the method, it was shown 
that, one could obtain the eigenf unctions and the eigenvalues of ES models simultaneously^^. 
The main aim of this present study is to demonstrate the applicability of the QHJ formalism 
to the periodic potentials. Since the QHJ formalism and its working are well studied in the 
literature, we do not go into details here but give only the necessary information. For details, 
the interested reader is referred to the earlier workJ^HSU an d the references there in. 

The main object of interest in the QHJ formalism is the quantum momentum function 



(QMF) p, which is the logarithmic derivative of the wavefunction ijj(x) } ' ' apart from a 
factor of —ih 

p = — ih— (In ip(x)). (5) 

(JjJb 

Hence, by obtaining the expression for the QMF one can obtain the expression for the 
wavefunction. The QMF satisfies the Riccati equation which is referred to as the QHJ 
equation, 

p 2 - ihp' = 2m(E - V(x)), (6) 
which is the special case of the general Riccati equation 

A(x)p 2 + B(x)p + C(x) + = 0. (7) 

Equation (JUJ) has the most convenient form to work with, in the QHJ formalism. Hence 
for all the cases studied through this method one tries to the bring the equation for the 
QMF to the form of (jOJ). In general, the solutions of Riccati equation (J2J) has two types 
of singularities, the fixed singularities and the moving singularities. The fixed singularities 
are determined by the singular points of A(x), B(x) and C(x). These appear in every 
solution and are independent of the initial conditions. Thus for the QHJ equation ©, the 
fixed singular points originate from the potential. The other singular points i. e, the moving 



singular points depend on the initial conditions and need not appear in every solution. It is 
known that for the Riccati equation only poles can appear as moving singularities ^1 

It is known that, the wavefunction for the n th excited state has n zeros. It is then seen 
from (JBJ that correspondingly p has n moving poles in between the two classical turning 
points and the residue at each of this pole is —%h. In general there could be other moving 
poles of p in the finite complex plane but one has little information regarding their location. 

In our earlier studies of the ES and QES solvable models it was assumed that the QMF 
has finite number of moving poles in the finite complex planep^Hl which turned out to be 
true for all the models studied. Therefore, for the study of the exactly solvable periodic 
potentials, without losing generality, one can assume that, the QMF has finite number of 
moving poles in the complex plane. It follows from our assumption that the point at infinity 
is an isolated singularity. 

With this assumption on the singularity structure of the QMF, we proceed to write the 
QHJ equation (JBJ) as 

q 2 + q' = V(x) - E, (8) 
with h = 2m = 1, where q = -^ip(x). The wavefunction in terms of q will be 

ijj(x) = exp( J q(x)dx). (9) 

One can see from (jHJ that the residue at the moving poles is unity. In general it has been 
found useful to change variable x to variable t = f(x) so as to make the coefficients of the 
Riccati equation rational. After a change of variable, (jHJ) becomes 

q 2 + F(t)^ + E-V(t) = 0, (10) 

where F(t) = 4- expressed as a function of t and V(t) is the potential in terms of t. One 
can see that the above equation does not have the convenient form of (JHJ). Hence to get (jlUJI 
into the form of (jHJ), we perform a transformation from q(x) to x(t) as follows, 

q = F(t)<P, ^ = x -ii(]nF(f)). (11) 

Using the above transformations (jlUj) becomes 

2 d X E-Vlt) l{F"(t)\ l{F'(t)\ 2 n 
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which is in the form of the Riccati equation. The residue at the moving poles is unity. The 
fixed poles correspond to the zeros of F(t). One can make use of (|T2*j) instead of the original 
QHJ equation for any general potential. In the next section we show, by taking the explicit 
example of the general Lame potential, how one can obtain the form of the band edge wave 
functions. 



III. GENERAL LAME POTENTIAL 



The QHJ equation in terms of q, for the Lame potential is 

q 2 + q' + E + l)msn 2 (x,m) = 0. (13) 

To proceed further we need to bring the potential to a meromorphic form, for which we do 
the following the change of variable 

t = sn (x, m), (14) 

with sn x = sn (x + 4K(m)). We would like to point out here that, the Lame equation can be 
written in five forms namely two algebraic, one trigonometric, one Weierstrassian and one 
Jacobff^, depending on the change of variable. Of all these change of variables we found 
()14|). which gives the Jacobi form of this Lame equation, to be best suited for our method. 
It enables one to write the QHJ equation in a form which can be easily analyzed using the 
QHJ formalism. Another added advantage of this particular choice is that, it maps half of 
the period parallelogram < x < 2K(m) of the Jacobi elliptic sn (x, m) function to the 
upper half complex planeJ^ One gets the equation for x as 

2 dx {mtf + 2m t 2 + 2 2E - 2j{j + l)mt 2 - mt 2 
X dt 4(1 -mt 2 ) 2 4(1 -t 2 ) 2 2(1 -t 2 ) (I -mt 2 ) ' 1 ' 



where 



i = V(i-* 2 )(i-™*% ^ = x+Kr^^ + rbO- (16) 

The following properties of the Jacobi elliptic functiond^l were used. 

-^-sn (x, m) = cn(x, m)dn(x, m) (17) 
dx 

and 

sn 2 (x, m) + cn 2 (x, m) = 1, sn 2 (x, m) + mdn 2 (x, m) = 1. (18) 



For all our further calculations we will regard x as the QMF instead of p. Using (|15jh 
in place of the original QHJ equation (j!3|) . one can obtain the expressions for the wave 
functions by analyzing the singularity structure of \ 

Singularity structure of \ : Equation (J15)) shows that \ has fixed poles at t — ±1 
and t = ±l/s/m and there are finite number of moving poles in the complex plane. Hence 
we make an assumption that the point at oo is an isolated singular point and that there are 
no singularities of the QMF except for those mentioned above. Therefore one can write Xi 
separating it into its singular and analytical parts in the following form, 

Here, Q(t) is the analytic part of x an d the rational terms represent the singular parts. 
Here, b^bi, d' x and d\ are the residues at t = ±1 and t = ±1/ ' \fm respectively. From ()15j) . 
one can see that x is bounded at infinity, which makes Q(t) in ()19p. analytic and bounded 
at infinity. Hence from Louville's theorem Q(t) will be a constant C. The summation term 
in (|19J). represents the singular part coming due to the finite number of moving poles for 
which the residues are easily found to be one and this term can be written as where 

* n 

n 

P n = \\{t — tk) is an n th degree polynomial. Thus (fT^j) becomes 



fc=0 



h b[ d x d[ P' n _ 

x = — + —i + T^iJ^ + T^iJ^i + T n +C - (20) 



The residues bx, b[, d\ and d\ at the fixed poles can be determined by substituting the 
Laurent expansion of x around these poles. For example, to calculate the residue at t — 1, 
one expands x 

X = ^Y + ao + a x (t-l) + (21) 

Substituting the above equation in (fTo]) and comparing the coefficients of different powers 
of (t — 1), one gets a quadratic equation in b\ whose roots are 

h = 3/4, 1/4. (22) 

Similarly the two values of the residues at all the other fixed poles b[, d\ and d[ turn out to 
be 

3/4, 1/4. (23) 
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Note that, all the values of the residues are independent of the potential parameter j. Hence, 
irrespective of the value of j in the potential, the residues at the fixed poles will take only 
the above values for any Lame potential. 

It will be shown later that, the knowledge of the residues is sufficient to obtain the wave 
function. Hence, at this stage, we should check, which values of the residues, out of the 
two values, give rise to acceptable' wave functions. In the present case, it turns out that 
both the values, i.e, 3/4 and 1/4 for each of the residues, leads to acceptable wave functions. 
The only requirement that restricts possible combinations of values is that of parity. The 
fact that if) has definite parity implies that x{~ t) = —x(t)i which rules out the possibilities 
bi 7^ b[ and d% ^ d[ and hence of all the possible combinations of the values of bi , b[ , d\ and 
d[, the only ones which are accepted are those with b\ = b' x and d\ = d' v 

In contrast to the above results, in our earlier study of ES boundstate problemJ^ an d 
QES problemJ^Sl one was lead to a unique choice of the residues at the fixed poles when one 
demanded that the wave function be square integrable. The other values of the residue did 
not correspond to any physical solution and were ruled out. In the case of potentials, which 
exhibit two phases of SUSY, for different ranges of the potential parameter both values 
of the residues were useful. One set of residues gave physically acceptable solutions for exact 
SUSY case and the other set gave for the broken SUSY. In the case of periodic potentials 
there is no such restriction as square integrability and so there is no way of ruling out one 
of the values. The reason that only one value of residue is acceptable for the bound state 
problems in one dimension can be attributed to the fact that the bound state solutions are 
non-degenerate. Whereas the band edge wavefunctions of the periodic potentials are doubly 
degenerate, we keep both the values of the residues. 

Behaviour of x a ^ infinity : We have assumed that the point at infinity is an isolated 
singularity. The form of (JT3J) suggests that x is bounded at oo. Therefore, x can be expressed 
as 



/ \ , Al Ao 

xW = A + - 1 + -f + .... 



(24) 



and the coefficients A^'s are fixed using the Riccati equation (JT5j) and one gets 



A, 



•o — 



(25) 



and two values for Ai as 



\x = 3 + 1, ~j- 



(26) 
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From (|26p. one obtains the leading behaviour of x{t) at oo which is dependent on j. Equation 
(J2UJ) gives the form of % which holds for all finite values of t. For large t, comparing (J2UJ) 
and (j21J), one gets 

2&i + 2di + n = A x . (27) 

As the left hand side is positive, it is clear that, only the choice Ai = j + 1 is consistent with 
the above equation, while the other value Ai = — j is ruled out. Thus from (|27jl . one gets, 
the expression for n as 

n = j + 1 - 2&i - 2di, (28) 

which gives the degree of the polynomial P n . 

Explicit forms of the wavefunctions : One can obtain the form of the wave functions 
using (JOJ). Substituting the relations given in (JTBJ) in (JOJ), one gets the expression for the 
wave function in terms of \ as ) 

*( ( ) = exp/( x ( i ) + i( r ^ + r l^)) rf( (29) 

rfil-Ab^t (l-4d 1 )mt P'\ , 

which when simplified and written in terms of the original variable x, gives 

ip(x) = (en x) a (da. XfP n (snx), (31) 

where, a = and /3 = . For the sake of simplicity, the elliptic modulus m of the 

Jacobi elliptic functions will be suppressed here onwards. The four different combinations 
of the residues (given as sets 1 to 4 ) give rise to four different forms of the band edge 
wavefunctions as listed in the table I. 

The parity constraint x{~ t) = —x(t), restricts the polynomial P n (t) to have either only 
odd or only even powers of t = snx. In the sixth and the seventh columns of table I, the 
number of linearly independent solutions, for the two cases, j being odd and j being even 
are given in terms of a positive integer N. For odd j, N = (j — l)/2 and for even j, N = j/2. 
In both the cases the total number of solutions for a particular j is equal to 2j + 1. It is 
easy to see that, the four forms and the number of solutions, of a particular form obtained 
here are in agreement with those already knowrP^I For a given set of &iand d\, n is fixed 
using (|28j) and the differential equation for the unknown polynomial p n (t) can be obtained 
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by substituting x{t) from (j2*U|l in (JT^J), which gives Q = and 

K + « + ^rjj PL + G{t)p n = o, (32) 

where 



G(t) 



t 2 {4b\ - 26i + 1/4) + 1/2 - 26i (mt) 2 (4d? - 2d x + 1/4) + m/2 - 2mdi 



(t 2 - 1) 



(mt 2 - l) 2 
2E - 2j(j + l)mt 2 + (16mMi - l)mt 2 
2(t 2 - l)(mt 2 - 1) 



For each set of residues given in table I, this differential equation is equivalent to a system of n 
linear equations, for the coefficients of different powers of t in P n (t). The energy eigenvalues 
are obtained by setting the corresponding determinant equal to zero. We illustrate this 
process by obtaining the eigenvalues and eigenfunctions explicitly for the case j = 2 



IV. LAME POTENTIAL WITH j = 2 

Using the procedure described in the previous section, we obtain the eigenvalues and the 
eigenfunctions for the supersymmetric potential 

VI (x) = 6m sn 2 x - 2m - 2 + 25, (33) 



where 5 = yl — m + m 2 . This potential is same as the Lame potential in (J2J) with j = 2, 
except for an additive constant 25 — 2m — 2, which has been added to make the lowest energy 
equal to zero. The equation for x{t) is, 

2 Y (mt) 2 + 2m t 2 + 2 2E + 4(m + 1 - 5) - 13mt 2 . , 

v — [— — - — [— — — ^ 1 1 i = 0. (34) 

X dt 4(1 -mt 2 ) 2 4(1 -t 2 ) 2 2(l-t 2 )(l-mt 2 ) 1 ' 

From (|28j) . the number of moving poles n, for j = 2 are 

n = 3 - 26i - 2di, (35) 

where the values of &i and di are obtained from (|2*2*|) and ()23)1 . For each set of b\, d\ and 
n, one can write the form and the number of linearly independent solutions by substituting 
j = 2, which gives iV = 1 in table I. Hence for the case j = 2, the form and the number 
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of solutions is as given in table II. To get the unknown polynomial part P n (t) of the wave 
function and the band edge energies, one needs to substitute the different sets of b x , d\ and 
n from table II, in the differential equation 

K + 4t + PL + G{t)P n = 0, (36) 

where 

_ ^ 2 ( 4fo ? - 2& i + 1/4) + 1/2 - 2&! (mt) 2 (4d 2 - 2d! + 1/4) + m/2 - 2mrfi 
W ~ (i 2 - l) 2 + {mt 2 - l) 2 + 

2£ + 4(m + 1 - 5) + (I6M1 - 13)mt 2 
2(t 2 -l)(mt 2 -l) ' 

Thus for various sets of residues one has the band edge energies and the wave functions 

as: 

Set 1 : b x = 1/4 , d x = 1/4, n = 2 
Taking P2 = At 2 + Bt + C, the parity constraint implies 

B = 0. (37) 

Substituting b 1; di and P 2 in (JSHJ), one gets a 2 x 2 matrix equation for A and C as follows 



0. (38) 



E- 2m -2 -25 -6m \ I A 

2 E + 2m + 2-25 ) ^0 

Equating the determinant of the matrix in (J38)) to zero, one gets the two values for energy 
as 

E I = AS, E u = 0, (39) 
which in turn give the polynomial as, 

Pj = m+ l-§- s m t 2 i Pjj = m + 1 + 5 - 3mt 2 . (40) 

From (|31jh we see that, the band edge wave functions in terms of x will be 

ipi(x) = m + 1 — 5 — 3m sn 2 x, ipu(x) = m + 1 + 5 — 3msn 2 (x). (41) 

Set 2 : b x = 3/4 , d x = 1/4, n = 1 
Taking Pi = t — tk and substituting in one gets the band edge eigenvalue and tk as, 

E = 25 -m + 2, t k = 0. (42) 
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The band edge wavef unction becomes 

ip(x) = cnisni. (43) 

Set 3 : b x = 1/4 , d x = 3/4, n = 1 
Proceeding in the same way as in set 2 one gets, 

E = 25 + 2m — 1, ip(x) = dnx snx. (44) 

Set 4 : 6i = 3/4 , di = 3/4, n = 
In this case Po as a constant. One gets the band edge energy and the wave function to be 

E = 25 — m — 1, ip(x) = cnx dnx. (45) 

Thus one obtains five band edge wave functions and their corresponding energies which agree 
with those given in^. 

V. ASSOCIATED LAME POTENTIAL 

The associated Lame potential is ES when a = b = j. In this case the potential expression 
becomes 

(cn ^ cc \ 
sn 2 xH . (46) 
dn 1/ 

The QHJ equation with the associated Lame potential is 

(cn 2 x \ 
sia 2 x H = °- ( 47 ) 
dn x J 

Doing the change of variable t = snx and proceeding in the same way as in the case of 
general Lame potential described in sec II, one gets the equation for x(t) as, 

2 dx (mt) 2 + 2m(l-2j(j + l)) 2 + t 2 2E - mt 2 (l + 2j(j + 1)) = 
X rft 4(1 -mt 2 ) 2 4(1 - 1 2 ) 2 2(1 -t 2 )(l -mt 2 ) ' 1 ] 

Singularity structure of x(t) : Similar to the Lame potential, x{t) nas fixed poles 
at t = ±1 and t = il/^/m, along with n moving poles in the entire complex t plane and 
further we assume that, the point at infinity is an isolated singularity, x has no other 
singular points, except at t = ±1 and t = ±l/-y/m. Following the same procedure, used in 
sec 77, one obtains the residues at t — ±1 as 

* = !•!■ « = !■!■ < 49 > 
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which is independent of j. The residues at t — ±.l/^fm turn out to be 

3 + 2j l-2j 3 + 2j l-2j 

which are j dependent. Knowing the singularity structure of x, one can write it as 

X t-l t+1 t-l/y/m t+l/y/m Pn 



valid for all t similar to the Lame potential. 

Behavior of x a t infinity : The behavior of x f° r large t, as fixed from the QHJ 
equation (jlHjl is 

X = A + y + ^ + .... (52) 
and one gets the two values for Ai as 

\i=j + 1, -j. (53) 

This result should agree with the large t behavior of x given by (|5T|). which is 

X = 6i + b[ + d x + d[ + n. (54) 

As before, parity requirement implies b\ = b^and d\ = d[ . Comparing the leading terms for 
large t, in (JoTj) and (|53j) we get) 

26i + 2di + n = A x . (55) 

In this case both the values of Ai are allowed because d± and d\ can take negative values 
unlike the case of Lame potential where \\ = —j was ruled out. Thus one has two cases 
Ai = j + 1 and Ax = — j. 
Case 1 : A = j + 1 

The values of n, which describe the number of zeros, for different combinations of b\ and 
d\ values are given in table III. The sets 3 and 4 give n = — 1 and n = —2 respectively and 
will not be considered, as n should be greater than or equal to 0. The sets 1 and 2 will give 
positive values of n only if j is positive. Hence, this table is used when j is positive. 

Case 2 : Ai = — j 

The n values for various sets of b\ and d\ with Ai = —j are listed in table IV . As in the 
previous case of Ai = j + 1, we will not consider the first two sets in the above table. 
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Note that the potential is invariant under the transformation j — > — j — 1. From the 
two tables III and IV, we see that two different values of j, one positive (j = j' > 0) 
and another negative (j = —j' — 1 < 0) leads to the same expression for the potential 
described by (@J). These two values also lead to the same answers for the wave functions and 
eigenvalues. Hence it is sufficient to restrict j to positive values alone. Thus using the two 
sets of combinations given in table III and the expressions for the wave functions given in 
(13 lj) one gets the following explicit forms and the number of solutions given in table V. 

In the next section, we obtain explicit expressions for the wavefunctions and energies for 
3=1- 



VI. ASSOCIATED LAME POTENTIAL WITH j = 1 

We perform the calculation with the supersymmetric potential 

cn 2 x 

VI (x) = 2m sn 2 x + 2m — * 2 — m + 2\/l — m, (56) 

dn x 

corresponding to j = 1. Proceeding in the same way as in the previous sections, the equation 
for x(t) is found to be 

2 dx (mtf-Qm t 2 + 2 2E + 4 + 2m - Ay/\ -m- 5mt 2 _ 
X + ~dt + A{mt 2 - I) 2 + 4(t 2 - l) 2 + 2(1 -mt 2 )(l - t 2 ) ~ ( ' 

and the form for x(t) is 

X t-1 t + 1 t-l/y/m t + l/y/m P n { ; 

Substituting (jH^j) in (|H7j) . one is left with the differential equation 

K + 4t (J^ + -S^rl K + G{t)P n = 0, (59) 



t 2 - 1 mt 2 - 1 



where 



_ t 2 (46 2 -2& 1 + l/4) + l/4-2fe 1 (mt) 2 (4rf 2 - 2di + 1/4) + -3m/2 - 2md 1 
(t} ~ (t 2 - l) 2 + (mt 2 - l) 2 + 

2E + 4 + 2m - Ay/1 - m + (166idi - 5)mt 



2 



2(i 2 - l)(mt 2 - 1) 

Substituting j = 1 in table III, one gets the band edge eigenfunctions and eigenvalues 
from set 1 and set 2 as follows. 
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Set 1 : 61 = 3/4 , di = -1/4, n = 1 

This combination gives only one solution with 



cn x sn x 



E = 2 - m + 2V1 - m, ib(x) = — : . (60) 

diax 



Set 2 : 61 = 1/4 , di = -1/4, n = 2 
This combination gives two solutions. The band edge energies are 



Ei = 0, E n = iVT^rn (61) 
and the corresponding wave functions are 



^/(x) = — dnx H , ^//(xj = — dnx , (62) 

m\ dnxy m \ dnx I 

which match with the wave functions i 



VII. CONCLUSIONS 



In this paper, we have studied the ES periodic potentials of the elliptic class namely the 
Lame and the associated Lame potentials. Using the QHJ formalism we have obtained here 
the general form of the band edge wave functions for these potentials, with the potential 
parameters taking only integer values . Also, the band edge eigenfunctions and the energy 
eigenvalues have been obtained explicitly for j = 2, for the Lame potential and p = q = 2, for 
the associated Lame potential. In the process, we have studied p, the logarithmic derivative 
of the wave function for these elliptic potentials and have found an interesting singularity 
structure exhibited by p in the complex domain. 

The QHJ formalism has been successful in obtaining the energy eigenvalues and the wave 
functions for a large variety of ES models and periodic potentials. The most important 
steps in obtaining the solutions using QHJ formalism have been, the choice of the change 
of variable and the ability to guess the right singularity structure of the QMF. It should be 
noted here that, our approach has been to assume a singularity structure of the QMF as 
simple as possible. Once a correct choice has been found, one could obtain the solutions in 
a most straight forward fashion even if the resulting equations were not of a simple form. 
From our earlier study of EsP*^, QES^l and the present study of periodic potentials, some 
interesting features of these models appear to be correct. All the models whether ES, QES 
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or periodic potentials, which we have studied so far, share a common property that 'the 
QMF', X-, becomes a rational function after a suitable change of variable. Therefore, in the 
finite complex plane it has finite number of moving poles, described by the parameter n. 
Our assumption that the point at infinity is an isolated singular point is equivalent to the 
condition that the QMF has finite number of moving poles. 

There are interesting differences in details of the singularity structure for different po- 
tential models. For the ES models, each value of n corresponds to an energy level and an 
eigenfunction with only n real zeros. These zeros correspond to the n moving poles of the 
QMF and are confined to the classical region only. For the QES models, n appears as a 
parameter in the potential and of the infinite possible states only n states can be obtained 
analytically. The QMF for all these states will have the same number of moving poles which 
are both real and complex. Correspondingly the wavefunctions for each state will have the 
same number (= n) of zeroes, of which the number of real zeros are in accordance with the 
oscillation theorem. The number of real zeros increases as one goes from the lowest state to 
the highest state possible. 

On the other hand, the exactly solvable periodic potential models show quite a different 
kind of distribution of the moving poles in the complex plane. Like QES models, the moving 
poles of the QMF consist of both real and complex poles. But unlike the QES models, where 
the number of poles of QMF for all the known states is same, for ES Periodic potential 
models one finds groups of solutions, with number of moving poles remaining same within a 
group but varying from one group to another group of solutions. This point becomes clear 
from the form of the solutions listed in table I for the Lame potential. Different groups of 
solutions are precisely the different sets of the solutions listed in table I. In each set, the 
degree of the polynomial P n (sn:r) is same and the total number of zeros of the band edge 
wave functions in the interval < x < 2K{m) are same. For example, for the N + 1 linear 
independent solutions belonging to set 1 the total number of zeros of the wavefunctions will 
be j. However for the solutions belonging to set 2 there will be j — 1 zeros of P n (sn x) and one 
zero corresponding to cn x, thus a total of j zeros. For all the sets the total number of zeros 
are given in the last column of table I. However the number of real zeros in a group increases 
with increasing energy and correspondingly the number of complex zeros will decrease. For 
all values of j, we will have utmost four such groups. A similar pattern is observed for ES 
associated Lame potential. 
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When one computes the residues at the fixed pole one gets two solutions due to the 
quadratic nature of the QHJ equation. For the bound state problems, both ES and QES 
considered in the earlier papersp^EH the square integrability of the wave function had to 
be insisted. This allowed us to pick the right residue, out of the two values, which gave the 
physically acceptable wave function, but in the present case of periodic potential models, 
there is no such restriction and, both values of the residues lead to acceptable solutions. The 
only restriction on the possible combination of the residues comes from the requirements of 
parity. 

To conclude we note that the new QHJ equation, obtained after a change of variable, 
involves rational functions in the independent variable. In order to make progress, we made 
crucial assumption that there are finite number of moving poles implying that the point at 
infinity is an isolated singular point. This important assumption has been used for all the 
models studied and its consistency with other equations has given useful restrictions on the 
acceptable solutions. Our present study completes this series of investigation of ES models 
from the QHJ approach. We now have a good starting point to take up models which are not 
ES or QES. For these models, making a proper guess about the singularity structure of the 
QMF does not appear to be possible. A study of the distribution of the moving poles in the 
complex plane may have a relation to the classical properties of the system. Such a study 
is interesting in its own right and can be best done numerically. We hope that here also, 
the established results in the complex variable theory will provide a useful scheme to obtain 
numerical, possibly approximate analytical solutions for the wave functions and energies. 
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TABLE I: Form of the wave functions i^(x), number of linear independent solutions and the total 
number of zeros of tp(x) for the general Lame potential V(x) = j(j + l)msn 2 x. 



sot 


bi 


di 


n 


in terms of x 


Number of solutions 


Total zeros 
of ^(x) 


j = 2iV + 1 


J = 2iV 


1 


1/4 


1/4 


3 


Pj(snx) 


N+ 1 


JV + 1 


i 


2 


3/4 


1/4 


i-i 


en x Pj- 1 (sni) 


7V+ 1 


A 


i 


3 


1/4 


3/4 


i-i 


dnxPj_i(snx) 


7V+ 1 


A 


i-i 


4 


3/4 


3/4 


i-2 


en xdn xPj-2 (sn x) 




A 


i-i 



TABLE II: The form of the wave functions and the number of linear independent solutions for the 
Lame potential with j = 2, V- = 6msn 2 x — 2m — 2 + 25 where S = y/1 — m + m 2 . 



set 


h 


di 


n 


ip in terms of x 


Number of LI solutions 


1 


1/4 


1/4 


2 


P 2 (snx) 


2 


2 


3/4 


1/4 


1 


cn xP\ (sn x) 


1 


3 


1/4 


3/4 


1 


dnxPi(snx) 


1 


4 


3/4 


3/4 





cn xdn xPq (sn x) 


1 



TABLE III: Values of n for different combinations of b\ and d\ for Ai = j + 1 for the associated 
Lame potential V(x) = j(j + l)m(sn 2 x + cn 2 x/dn 2 x). 



set 


6i 


di 


n 


1 


1/4 


(1 " 2i)/4 


2i 


2 


3/4 


(1 " 2j)/4 


2J-1 


3 


1/4 


(3 + 2j)/4 


-1 


4 


3/4 


(3 + 2j)/4 


-2 
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TABLE IV: Various sets of 61 and d\ for Ai = —j for associated Lame potential V(x) = j(j + 
l)m(sn 2 x + cn 2 x/dn 2 x). 



set 




di 


n 


1 


1/4 


(1 " 2i)/4 


-1 


2 


3/4 


(1 " 2i)/4 


-2 


3 


1/4 


(3 + 2j)/4 


-2j - 2 


4 


3/4 


(3 + 2j)/4 


-2j - 3 



TABLE V: Form of the wavefunctions for the associated Lame potential V(x) = j(j + l)m(sn 2 x + 
cn 2 x/dn 2 x) for positive values of j. 



set 


form of ip 


No of solutions 






j = 2N 


j = 2N + 1 


1 


cnxP2j-i(snx) 


2N 


2N + 1 


(dnx)i 


2 


P 2 j(snx) 
(dnx)i 


2N + 1 


2N + 2 
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